arXiv: 1504.02025vl [hep-ex] 8 Apr 2015 


Measurements of ^’( 3686 ) —>■ K AH+ + c.c. and - 1 /^( 3686 ) —>• 'yK AH+ + c.c. 


M. Ablikim\ M. N. Achasov®’“, X. C. AiS O. Albayrak®, M. Albrecht'^, D. J. Ambrose'^^, A. Amoroso'^®^'^®*^', F. F. AnS 
Q. An**®, J. Z. Bai^ R. Baldini Ferroli^®'^, Y. Ban®\ D. W. Bennett^®, J. V. Bennett®, M. Bertani2°^, D. Bettoni^^'^, 

J. M. Bian^®, F. Bianchi^®'^’'‘®‘®', E. Boger^®’®^, O. Bondarenko^®, L Boyko^®, R. A. Briere®, H. Cai®®, X. Cai®, O. Cakir'^®'^’®, 
A. Calcaterra^®^, G. F. Cao®, S. A. Cetin^°®, J. F. Chang\ G. Ghelkov^®’^ G. Ghen®, H. S. Chen\ H. Y. Ghen^, 

J. G. Chen\ M. L. Chen\ S. J. Chen^®, X. Chen\ X. R. Chen^®, Y. B. Ghen®, H. P. Gheng®’’, X. K. Ghu®®, G. Cibinetto^®'^, 
D. Gronin-Hennessy^®, H. L. Dai^, J. P. Dai®"*, A. Dbeyssi^^, D. Dedovich^®, Z. Y. Deng^, A. Denig^^, I. Denysenko^®, 

M. Destefanis"'®'^-'*®'^, F. De Mori'‘®'^’‘‘®'^, Y. Ding^'^, C. Dong®°, J. Dong\ L. Y. Dong®, M. Y. Dong®, S. X. Du®^ 

P. F. Duan®, J. Z. Fan®®, J. Fang®, S. S. Fang®, X. Fang"®®, Y. Fang®, L. Fava®®®’®®®®, F. Feldbauer®®, G. Felici®®"®, 

G. Q. Feng®®, E. Fioravanti^®^, M. Fritsch®®’^^ G. D. Fu®, Q. Gao®, X. Y. Gao®, Y. Gao®®, Z. Gao®®, I. Garzia®®^, 

C. Geng®®, K. Goetzen®°, W. X. Gong®, W. Gradl®®, M. Greco®®®®’®®*^, M. H. Gu®, Y. T. Gu®®, Y. H. Guan®, A. Q. Guo®, 

L. B. Guo®®, Y. Guo®, Y. P. Guo®®, Z. Haddadi®®, A. Hafner®®, S. Han®°, Y. L. Han®, X. Q. Hao®®, F. A. Harris®®, K. L. He®, 

Z. Y. He®°, T. Held®, Y. K. Heng®, Z. L. Hou®, G. Hu®®, H. M. Hu®, J. F. Hu®®®®’®®'^, T. Hu®, Y. Hu®, G. M. Huang®, 

G. S. Huang®®, H. P. Huang®®, J. S. Huang®®, X. T. Huang®®, Y. Huang®®, T. Hussain®'^, Q. Ji®, Q. P. Ji®°, X. B. Ji®, 

X. L. Ji®, L. L. Jiang®, L. W. Jiang®®, X. S. Jiang®, J. B. Jiao®®, Z. Jiao®^, D. P. Jin®, S. Jin®, T. Johansson®®, A. Julin®®, 

N. Kalantar-Nayestanaki®®, X. L. Kang®, X. S. Kang®®, M. Kavatsyuk®®, B. C. Ke®, R. Kliemt®®, B. Kloss®®, 

O. B. Kolcu®®®’'', B. Kopf®, M. Kornicer®®, W. Kiihn®®, A. Kupsc®®, W. Lai®, J. S. Lange®®, M. Lara®®, P. Larin®®, 

G. Leng®®®^, C. H. Li®, Cheng Li®®, D. M. Li®®, F. Li®, G. Li®, H. B. Li®, J. G. Li®, Jin Li®®, K. Li®®, K. Li®®, Lei Li®, 

P. R. Li®®, T. Li®®, W. D. Li®, W. G. Li®, X. L. Li®®, X. M. Li®®, X. N. Li®, X. Q. LP, Z. B. Li®®, H. Liang®®, Y. F. Liang®®, 

Y. T. Liang®®, G. R. Liao®®, D. X. Lin®®, B. J. Liu®, G. X. Liu®, F. H. Liu®®, Fang Liu®, Feng Liu®, H. B. Liu®®, H. H. Liu®, 

H. H. Liu®®, H. M. Liu®, J. Liu®, J. P. Liu®°, J. Y. Liu®, K. Liu®^, K. Y. Liu®'', L. D. Liu®®, P. L. Liu®, Q. Liu®®, S. B. Liu®®, 

X. Liu®®, X. X. Liu®®, Y. B. Liu®°, Z. A. Liu®, Zhiqiang Liu®, Zhiqing Liu®®, H. Loehner®®, X. G. Lou®’% H. J. Lu®', 

J. G. Lu®, R. Q. Lu®®, Y. Lu®, Y. P. Lu®, G. L. Luo®®, M. X. Luo®®, T. Luo®®, X. L. Luo®, M. Lv®, X. R. Lyu®®, F. G. Ma®', 
H. L. Ma®, L. L. Ma®®, Q. M. Ma®, S. Ma®, T. Ma®, X. N. Ma®°, X. Y. Ma®, F. E. Maas®®, M. Maggiora®®"®’®®'®', 

Q. A. Malik®’', Y. J. Mao®®, Z. P. Mao®, S. Marcello®®'®’®®'^, J. G. Messchendorp®®, J. Min®, T. J. Min®, R. E. Mitchell®®’, 
X. H. Mo®, Y. J. Mo®, C. Morales Morales®®, K. Moriya®®*, N. Yu. Muchnoi®’’“, H. Muramatsu®®, Y. Nefedov®®, F. Nerling®®, 

I. B. Nikolaev®’’“, Z. Ning®, S. Nisar®, S. L. Niu®, X. Y. Niu®, S. L. Olsen®®, Q. Ouyang®, S. Pacetti®°®, P. Patteri®®^, 

M. Pelizaeus®, H. P. Peng®®, K. Peters®®, J. Pettersson®®, J. L. Ping®®, R. G. Ping®, R. Poling®®, Y. N. Pu®®, M. Qi®®, 

S. Qian®, C. F. Qiao®®, L. Q. Qin®®, N. Qin®®, X. S. Qin®, Y. Qin®®, Z. H. Qin®, J. F. Qiu®, K. H. Rashid®', C. F. Redmer®®, 

H. L. Ren®®, M. Ripka®®, G. Rong®, X. D. Ruan®®, V. Santoro®®'®, A. Sarantsev®®’-^, M. Savrie®®®, K. Schoenning®®, 

S. Schumann®®, W. Shan®®, M. Shao®®, G. P. Shen®, P. X. Shen®°, X. Y. Shen®, H. Y. Sheng®, W. M. Song®, X. Y. Song®, 

S. Sosio®®^-®®'®', S. Spataro®®'®'®®'®', G. X. Sun®, J. F. Sun®®, S. S. Sun®, Y. J. Sun®®, Y. Z. Sun®, Z. J. Sun®, Z. T. Sun®®, 

G. J. Tang®®, X. Tang®, 1. Tapan®®®^, E. H. Thorndike®®, M. Tiemens^®, D. Toth®®, M. Ullrich^®, 1. Uman®°®, G. S. Varner®®, 
B. Wang®°, B. L. Wang®®, D. Wang®®, D. Y. Wang®®, K. Wang®, L. L. Wang®, L. S. Wang®, M. Wang®®, P. Wang®, 

P. L. Wang®, Q. J. Wang®, S. G. Wang®®, W. Wang®, X. F. Wang®®, Y. D. Wang®°^, Y. F. Wang®, Y. Q. Wang®®, Z. Wang®, 

Z. G. Wang®, Z. H. Wang®®, Z. Y. Wang®, T. Weber®®, D. H. Wei®®, J. B. Wei®®, P. Weidenkaff®®, S. P. Wen®, U. Wiedner®, 
M. Wolke®®, L. H. Wu®, Z. Wu®, L. G. Xia®®, Y. Xia®«, D. Xiao®, Z. J. Xiao®®, Y. G. Xie®, Q. L. Xiu®, G. F. Xu®, L. Xu®, 
Q. J. Xu®®, Q. N. Xu®®, X. P. Xu®', L. Yan®®, W. B. Yan®®, W. C. Yan®®, Y. H. Yan®®, H. X. Yang®, L. Yang®°, Y. Yang®, 

Y. X. Yang®®, H. Ye®, M. Ye®, M. H. Ye', J. H. Yin®, B. X. Yu®, G. X. Yu®°, H. W. Yu®®, J. S. Yu^®, G. Z. Yuan®, 

W. L. YuaP®, Y. Yuan®, A. Yuncu®°®’®, A. A. Zafar®', A. Zallo^®'®, Y. Zeng®®, B. X. Zhang®, B. Y. Zhang®, G. Zhang^®, 

G. C. Zhang®, D. H. Zhang®, H. H. Zhang®®, H. Y. Zhang®, J. J. Zhang®, J. L. Zhang®, J. Q. Zhang®, J. W. Zhang®, 

J. Y. Zhang®, J. Z. Zhang®, K. Zhang®, L. Zhang®, S. H. Zhang®, X. Y. Zhang®®, Y. Zhang®, Y. H. Zhang®, Y. T. Zhang®®, 

Z. H. Zhang®, Z. P. Zhang®®, Z. Y. Zhang®®, G. Zhao®, J. W. Zhao®, J. Y. Zhao®, J. Z. Zhao®, Lei Zhao®®, Ling Zhao®, 

M. G. Zhao®°, Q. Zhao®, Q. W. Zhao®, S. J. Zhao®^ T. C. Zhao®, Y. B. Zhao®, Z. G. Zhao®®, A. Zhemchugov^®’'*, B. Zheng®®, 

J. P. Zheng®, W. J. Zheng®®, Y. H. Zheng®®, B. Zhong^®, L. Zhou®, Li Zhou®°, X. Zhou®®, X. K. Zhou®®, X. R. Zhou®®, 

X. Y. Zhou®, K. Zhu®, K. J. Zhu®, S. Zhu®, X. L. Zhu®®, Y. C. Zhu®®, Y. S. Zhu®, Z. A. Zhu®, J. Zhuang®, L. Zotti®®^’®®®', 

B. S. Zou®, J. H. Zou® 


(BESIII Collaboration) 

® Institute of High Energy Physics, Beijing 10004-9, People’s Republic of China 
® Beihang University, Beijing 100191, People’s Republic of China 
® Beijing Institute of Petrochemical Technology, Beijing 10261 7, People’s Republic of China 
® Bochum Ruhr-University, D-jfISO Bochum, Germany 
® Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA 
® Central China Normal University, Wuhan 430079, People’s Republic of China 
' China Center of Advanced Science and Technology, Beijing 100190, People’s Republic of China 
COMSATS Institute of Information Technology, Lahore, Defence Road, Off Raiwind Road, 54000 Lahore, Pakistan 
® G.I. Budker Institute of Nuclear Physics SB RAS (BINP), Novosibirsk 630090, Russia 
®® GSI Helmholtzcentre for Heavy Ion Research GmbH, D-64291 Darmstadt, Germany 
®® Guangxi Normal University, Guilin 541004, People’s Republic of China 
®® GuangXi University, Nanning 530004, People’s Republic of Ghina 
®® Hangzhou Normal University, Hangzhou 310036, People’s Republic of China 


2 


Helmholtz Institute Mainz, Johann-Joachim-Becher-Weg 45, D-55099 Mainz, Germany 
Henan Normal University, Xinxiang 453007, People’s Republic of China 
Henan University of Science and Technology, Luoyang 471003, People’s Republic of China 
Huangshan College, Huangshan 245000, People’s Republic of China 
Hunan University, Changsha 410082, People’s Republic of China 
Indiana University, Bloomington, Indiana 47405, USA 

(A)INFN Laboratori Nazionali di Frascati, 1-00044! Frascati, Italy; (B)INFN and University of Perugia, 1-06100, Perugia, 

Italy 

(A)INFN Sezione di Ferrara, 1-44132, Ferrara, Italy; (B)University of Ferrara, 1-44132, Ferrara, Italy 
Johannes Gutenberg University of Mainz, Johann-Joachim-Becher-Weg 45, D-55099 Mainz, Germany 
Joint Institute for Nuclear Research, 141980 Dubna, Moscow region, Russia 
Justus Liebig University Giessen, II. Physikalisches Institut, Heinrich-Buff-Ring 16, D-35392 Giessen, Germany 
KVI-CART, University of Groningen, NL-9747 AA Groningen, The Netherlands 
Lanzhou University, Lanzhou 730000, People’s Republic of China 
Liaoning University, Shenyang 110036, People’s Republic of China 
Nanjing Normal University, Nanjing 210023, People’s Republic of China 
Nanjing University, Nanjing 210093, People’s Republic of China 
Nankai University, Tianjin 300071, People’s Republic of China 
Peking University, Beijing 100871, People’s Republic of China 
Seoul National University, Seoul, 151-747 Korea 
Shandong University, Jinan 250100, People’s Republic of China 
Shanghai Jiao Tong University, Shanghai 200240, People’s Republic of China 
Shanxi University, Taiyuan 030006, People’s Republic of China 
Sichuan University, Chengdu 610064, People’s Republic of China 
Soochow University, Suzhou 215006, People’s Republic of China 
Sun Yat-Sen University, Guangzhou 510275, People’s Republic of China 
Tsinghua University, Beijing 100084, People’s Republic of China 
(A) Istanbul Ay din University, 34395 Sefakoy, Istanbul, Turkey; (B)Dogus University, 34722 Istanbul, Turkey; (C)Uludag 

University, 16059 Bursa, Turkey 

University of Chinese Academy of Sciences, Beijing 100049, People’s Republic of China 
University of Hawaii, Honolulu, Hawaii 96822, USA 
University of Minnesota, Minneapolis, Minnesota 55455, USA 
University of Rochester, Rochester, New York 14627, USA 
University of Science and Technology of China, Hefei 230026, People’s Republic of China 
University of South China, Hengyang 421001, People’s Republic of China 
University of the Punjab, Lahore-54590, Pakistan 

(A)University of Turin, 1-10125, Turin, Italy; (B)University of Eastern Piedmont, 1-15121, Alessandria, Italy; (C)INFN, 

1-10125, Turin, Italy 

Uppsala University, Box 516, SE-75120 Uppsala, Sweden 
Wuhan University, Wuhan 430072, People’s Republic of China 
Zhejiang University, Hangzhou 310027, People’s Republic of China 
Zhengzhou University, Zhengzhou 450001, People’s Republic of China 

“ Also at the Novosibirsk State University, Novosibirsk, 630090, Russia 
® Also at Ankara University, 06100 Tandogan, Ankara, Turkey 
“ Also at the Moscow Institute of Physics and Technology, Moscow 141700, Russia and at the Functional Electronics 

Laboratory, Tomsk State University, Tomsk, 634050, Russia 
Currently at Istanbul Arel University, 34295 Istanbul, Turkey 
® Also at University of Texas at Dallas, Richardson, Texas 75083, USA 
^ Also at the NRC ’’Kurchatov Institute”, PNPI, 188300, Gatchina, Russia 
® Also at Bogazici University, 34342 Istanbul, Turkey 
^ Also at the Moscow Institute of Physics and Technology, Moscow 14 1 700, Russia 


1 


Using a sample of 1.06 x 10® ip{S686) events produced in e''"e“ collisions at ^/s — 3.686 GeV 
and collected with the BESIII detector at the BEPCII collider, we present studies of the decays 
i/>(3686) —>■ K~A3'^ + c.c. and ip{36S6) —>■ •yK~NE'^ + c.c.. We observe two hyperons, H(1690)“ 
and H(1820)“, in the K~A invariant mass distribution in the decay ip{3686) —>■ K~AE'^ + c.c. with 
significances of 4.9a and 6.2a, respectively. The branching fractions of ip{3686) —>■ K~AE'^ + c.c., 
V’(3686) —>■ + c.c., i/>(3686) 7 XcJ 'yK~AE'^ + c.c. (J = 0, 1, 2), and ip{3686) —>■ 



3 


H(1690/1820) + c.c with subsequent decay H(1690/1820) K A are measured for the first 

time. 

PACS numbers: 13.25.Gv, 13.30. Eg, 14.20.Jn 


I. INTRODUCTION 

The quark model, an outstanding achievement of the 
last century, provides a rather good description of the 
hadron spectrum. However, baryon spectroscopy is far 
from complete, since many of the states expected in the 
SU(3) multiplets are either undiscovered or not well es¬ 
tablished [1], especially in the case of cascade hyperons 
with strangeness S = —2, the S*. Due to the small pro¬ 
duction cross sections and the complicated topology of 
the final states, only eleven S* states have been observed 
to date. Few of them are well established with spin-parity 
determined, and most observations and measurements to 
date are from bubble chamber experiments or diffractive 
K~p interactions [2]. 

As shown by the Particle Data Group (PDG), most 
S* hyperon results obtained to date have limited statis¬ 
tics [2]. For example, the S(1690)“ was first observed in 
the T,K final state in the reaction K~p —>■ {YjK)Ktt [3]. 
Afterwards its existence has been confirmed by other ex¬ 
periments [4-6], but its spin-parity was not well deter¬ 
mined. More recently, BABAR reported evidence for 
= 1/2“ for the S(1690) by analyzing the Legendre 
Polynomial moments of the S“7r“'" system in the decay 
A/" —>■ [7]. Clear evidence for S(1820) was 

observed in the AT “A mass spectrum from a sample of 
130±16 events in K~p interactions [8], and the J = 1/2 
assumption was ruled out by using the Byers and Fenster 
technique [9]. Ten years later, a CERN-SPS experiment 
indicated that S(1820) favors negative parity in the case 
of J = 3/2 [10]. 

At present, the S(1690) and S(1820) are firmly es¬ 
tablished. Further investigation of their properties, e.g. 
mass, width and spin-parity, is important to the under¬ 
standing of S* states. Besides scattering experiments, 
decays from charmonium states offer a good opportunity 
to search for additional S* states. Although charmonium 
decays into pairs of states are suppressed by the lim¬ 
ited phase space, the narrow charmonium width which 
reduces the overlap with the neighboring states and the 
low background allow the investigation of these hyperons 
with high statistics charmonium samples. 

Furthermore, our knowledge of charmonium decays 
into hadrons, especially to hyperons, is limited. The pre¬ 
cise measurements of the branching fractions of charmo¬ 
nium decays may help provide a better understanding 
of the decay mechanism. The large 'tp{3686) data sam¬ 
ple collected with the BESIII detector provides a good 
opportunity to study the cascade hyperons. 

In this paper, we report on a study of the decays 
ip{3686) — >■ K~AE^ +C.C. and'i/'(3686) —>■ 'yK~AE'^ + c.c. 
based on a sample of 1.06x10® ■i/'(3686) events [11] col¬ 
lected with the BESIII detector. Another data sample. 


consisting of an integrated luminosity of 44.5 pb“^ [12] 
taken below the tp{3686) peak at ^/s = 3.65 GeV, is used 
to estimate continuum background. Evidence for the 
S(1690)“ and S(1820)“ is observed in the K~A invariant 
mass distribution in the decay ■i/'(3686) —>■ AT” AS+ -|- c.c. 
In the following, the charge conjugate decay mode is al¬ 
ways implied unless otherwise specified. 


II. DETECTOR AND MONTE CARLO 
SIMULATION 

BEPCII is a two-ring collider designed for a luminosity 
of 10®® cm“^s“® at the ip{3770) resonance with a beam 
current of 0.93 A. The BESIII detector has a geometri¬ 
cal acceptance of 93% of 47r, and consists of a helium-gas- 
based drift chamber (MDG), a plastic scintillator time-of- 
flight system (TOF), a CsI(Tl) electromagnetic calorime¬ 
ter (EMC), a superconducting solenoid magnet provid¬ 
ing 1.0 T magnetic field, and a resistive plate chamber- 
based muon chamber (MUG). The momentum resolution 
of charged particles at 1 GeV/c is 0.5%. The time resolu¬ 
tion of the TOF is 80 ps in the barrel detector and 110 ps 
in the end cap detectors. The photon energy resolution at 
1 GeV is 2.5% (5%) in the barrel (end caps) of the EMC. 
The trigger system is designed to accommodate data tak¬ 
ing at high luminosity. A comprehensive description of 
the BEPCII collider and the BESIII detector is given in 
Ref. [13]. 

A GEANT4-based [14] MG simulation software 
BOOST [15], which includes geometric and material de¬ 
scription of the BESIII detector, detector response and 
digitization models as well as tracking of the detector 
running condition and performance, is used to gener¬ 
ate MC samples. A series of exclusive MG samples, 
V’(3686) ^ lXcJ_^ xK-AE+, V'(3686) ^ R-AE+, 
1/(3686) —>■ A'“E®S+ are generated to optimize the selec¬ 
tion criteria and estimate the corresponding selection ef¬ 
ficiencies. The production of i/(3686) is simulated by the 
generator KKMC [16, 17]. The decay i/(3686) —>■ xXcJ 
is assumed to be a pure El transition and to follow a 
1 + acos^O angular distribution with a = 1,—1/3 and 
1/13 for J = 0, 1 and 2, respectively [18], where 9 is the 
polar angle of the photon. The other subsequent decays 
are generated with BesEvtGen [19] with a uniform distri¬ 
bution in phase space. An inclusive MC sample, consist¬ 
ing of 1.06X 10® ■0(3686) events, is used to study potential 
backgrounds, where the known decay modes of 0(3686) 
are generated by BesEvtGen with branching fractions at 
world average values [2], and the remaining unknown de¬ 
cay modes are modeled by LUNDCHARM [20]. 
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III. ANALYSIS OF V’(3686) ^ R-AS+ 

The decay ip{3686) —>■ A“AS+ is reconstructed from 
the cascade decays A — >• ptt~, S + —>• A7r+ and A —>• p7r+. 
At least six charged tracks are required and their polar 
angles 9 must satisfy | cosd| < 0.93. The combined TOF 
and dE/dx information is used to form particle identifi¬ 
cation (PID) confidence levels for pion, kaon and proton 
hypotheses. Each track is assigned to the particle hy¬ 
pothesis type with the highest confidence level. Candi¬ 
date events are required to have one kaon. If more than 
one kaon candidate is identified, only the kaon with high¬ 
est confidence level is kept, and the others are assumed 
to be pious. The same treatment is implemented for the 
proton and antiproton. The final identified charged kaon 
is further required to originate from the interaction point 
(IP), i.e., the point of its closest approach to the beam 
is within 1 cm in the plane perpendicular to beam and 
within ±10 cm along the beam direction. 

In the analysis, constraints on the secondary decay ver¬ 
tices of the long lived particles, A and are utilized 
to suppress backgrounds. A particles are reconstructed 
using secondary vertex fits on pTr~ pairs. For events with 
more than one A candidate, the one with the smallest 
for the secondary vertex fit is selected. candidates 
are reconstructed in two steps. A p7r+ pair sharing a 
common vertex is selected to reconstruct the A candi¬ 
date, and the common vertex is regarded as its decay 
vertex. The is then reconstructed with a A candidate 
and another 7r+ by implementing another secondary ver¬ 
tex fit. For events with more than one candidate, the 
p7r+7r“'' combination with the minimum |M(p7r+) —M(A)| 
is selected, where M(p7r+) is the invariant mass of the A 
candidate from the secondary vertex fit, and M{A) is the 
corresponding nominal mass from the PDG [2]. 

The selected K~, A, and candidates are subjected 
to a four-momentum constraint kinematic fit (4C-fit) un¬ 
der the hypothesis of ■i/'(3686) —>■ K~ hEA, and x\c ^ 
is required to further suppress the potential backgrounds 
and to improve the resolution. Figure 1 (a) shows the 
invariant mass distribution of P'k~ , M{pTr~), where a 
A peak is clearly visible. A mass window requirement 
1.110 < M{pTr~) < 1.121 GeV/c^, corresponding to 6 
times the mass resolution, is imposed to select A can¬ 
didates. With the above selection criteria, the invari¬ 
ant mass of the A candidate M(p7r+) is shown in Fig. 1 
(b), and a clean A peak is observed. A mass window 
requirement 1.110 < M(p7r+) < 1.121 GeV/c^ is ap¬ 
plied to further improve the purity. Figure 1 (c) shows 
the scatter plot of M{jj'k~ ) versus without the 

A mass window requirement, where the accumulated 
events around the A-5 mass region are from the decay 
^(3686) — >■ K~KEA. The projection of M(A7r+) for all 
surviving events is shown in Fig. 1 (d), where the S’*' 
peak is seen with very low background. 

Potential non-5+ backgrounds are studied with the 
^/’(3686) inclusive MC sample by imposing the same selec¬ 
tion criteria. The corresponding distribution of M(A7r+) 


is shown in Fig. 1 (d) as the shaded histogram. The 
background is well described by the inclusive MG sample 
and is flat. Backgrounds are also investigated with the 
M{pTT~) versus M(p7r+) 2-dimensional sideband events 
from the data sample, where the sideband regions are 
defined as 1.102 < M(p7r“/p7r+) < 1.107 GeV/c^ and 
1.124 < M{p-K~ /pTT^) < 1.130 GeV/c^. No peaking 
structure is observed in the M{A'k~^) distribution around 
the region. To estimate the non-resonant background 
coming directly from e+e“ annihilation, the same selec¬ 
tion criteria are implemented on the data sample taken 
at ^/s = 3.65 GeV. Only 1 event with Af(A7r+) at 1.98 
GeV/c^, located outside of the signal region, survives, 
which is normalized to an expectation of 3.6 events in 
^p(S686) data after considering the integrated luminosi¬ 
ties and an assumed 1 /s dependence of the cross section, 
as L('s/s) oc A^obs/o’QED(-\/s) [ll]j where L is the inte¬ 
grated luminosity and ctqed is the cross section of QED 
processes. Therefore, the non-resonant background can 
be neglected. 

A. BRANCHING FRACTION MEASUREMENT 

To determine the event yield, an extended unbinned 
maximum likelihood fit is performed on the M (Att’’') dis¬ 
tribution in Fig. 1 (d). In the fit, the S"*" is described by a 
double Gaussian function, and the background is param¬ 
eterized by a first order Ghebychev polynomial function. 
The fit result, shown as the solid curve in Fig. 1 (d), 
yields Aobs = 236.4 ± 16.6 candidates. The decay 
branching fraction B{tj;{3686) K~AEA) is calculated 
to be 

^(V'(3686) ^ A:-AS+) 

_ _ Nohs _ 

-^i/>( 3686 ) ■ B‘^{A —>■ pTT~) ■ —>■ Att”) • e 

= (3.86 ±0.27) X 10-^ (1) 

where A^^( 3686 ) = (106.41 ± 0.86) x 10® is the number 
of i/'( 3686) events determined with inclusive hadronic 
events [11], e = 14.1% is the detection efficiency, evalu¬ 
ated from the MG sample simulated with a uniform dis¬ 
tribution in phase-space, and B{A —>■ pTr~) and B{E~^ —>■ 
Att’’') are the corresponding decay branching fractions [2]. 
The uncertainty is statistical only. 

B. OBSERVATION OF H*" STATES 

In the distribution of the K~A invariant mass, 
M{K~A), structures around 1690 and 1820 MeV/c^, as¬ 
sumed to be S(1690)” and S(1820)“, are evident with 
rather limited statistics. In order to improve the statis¬ 
tics, a partial reconstruction method is used where the 
K~ and A are required but the reconstruction of 
and the 4G kinematic fit are omitted. In addition, an 
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FIG. 1. Invariant mass distribution of (a) p7r“ and (b) (with the A mass window requirement). The arrows indicate the 
mass windows used in the analysis (see text), (c) Scatter plot of M{p'K~) versus M{A'k'^) for data, (d) Att"*" invariant mass 
distribution. In the one dimensional plots, the points with error bars are the data, the solid histograms are MC distributions 
normalized to the data, and the shaded histogram is the background estimated from the inclusive MC sample. The solid and 
long-dashed lines represent the fit curve and the background contribution from the fit. 


identified anti-proton is required among the remaining 
charged tracks to suppress background. With the above 
loose selection criteria, the distribution of M{pTr~) is 
shown in Fig. 2 (a), where a A is observed. After ap¬ 
plying the M{pTT~) mass window requirement, 1.110 < 
M{p'k~) < 1.121 GeV/c^, the distribution of the mass 
recoiling against the K~K system RM{K~A) is shown 
in Fig. 2 (b), where the is observed, although with a 
higher background than in the full reconstruction. With 
a requirement of 1.290 < RM{K~A) < 1.345 GeV/c^, 
theS(1690)“ andS(1820)“ are observed in the M{K~ A) 
distribution with improved statistics, as shown in Fig. 3. 
MC studies show that the event selection efficiency is im¬ 
proved by a factor of two using the partial reconstruction 
method. 

To ensure that the observed structures are not from 
background, potential backgrounds are investigated us¬ 
ing both data and inclusive MC samples. Non-A (S+) 
background is estimated from the events in the A (S+) 
sideband regions, defined as 1.102 < M{pTr~) < 1.107 
GeV/c2 and 1.124 < M{ptt-) < 1.130 GeV/c^ (1.243 < 


RM{K-A) < 1.270 GeV/c^ and 1.365 < RM{K-A) < 
1.393 GeV/c^), and their M{K~K) distribution is shown 
in Fig. 3 with the dot-dashed (dashed) histogram. Pos¬ 
sible background sources are also investigated with the 
inclusive MC sample, and the result is shown with the 
shaded histogram in Fig. 3. No evidence of peaking struc¬ 
tures in the M{K~K) distribution is observed in either 
the sideband region or the inclusive MC sample. The 
same selection criteria are applied to the data sample col¬ 
lected at 3.65 GeV to estimate the background coming 
directly from e“''e“ annihilation. Only one event with 
M{K~A) around 1.98GeV survives, which corresponds 
to an expected 3.6 events when normalized to the ■0(3686) 
sample. This background can therefore be neglected. 

An extended unbinned maximum likelihood fit of the 
M{K~ A) distribution is performed to determine the res¬ 
onance parameters and event yields of the excited 
hyperons. In the ht, the shapes are described by 
Breit-Wigner functions Ai (to) convoluted with Gaussian 
functions G{m, p,a), which represent the mass shift and 
resolution in the reconstruction, multiplied by the mass 
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RM(K A) (GeV/c^) 


FIG. 2. Invariant mass spectrum (a) of p 7 r“, and (b) of the mass recoiling against the K~A system. The dots with error bars 
show the distribution for data, and the solid histogram shows that for the exclusive MC normalized to the data in the signal 
region. The arrows indicate the selection region used in the analysis (see text). 



M(K'A) (GeV/c^) 

FIG. 3. Invariant K~A mass distribution. Points with er¬ 
ror bars represent data, and the solid and dashed curves are 
the fit curve and the non-resonant contribution obtained from 
the fit. The shaded histogram represents the background es¬ 
timated from the inclusive MC sample, and the dashed and 
dot-dashed histograms are the A sideband and the S'*" side¬ 
band backgrounds from data, respectively. 


dependent efficiency s{m), e(m)[G{m, 0 Af (m)]. 

In the fit, both parameters of and e(m) are 

fixed the the values determined from the studies to exclu¬ 
sive MC samples, and the Breit-Wigner function Ai{m) 
is described below. The shape of background is parame¬ 
terized by a function i?(m) = (m— 77 io)^^^-|-c(m— 
where mg is the mass threshold and c is a free parameter. 

The Breit-Wigner function Ai{m) used in the fit can 
be written as 

A(m) = -==-=- • 

m — M + 

I r 

where M, T are the mass and width of the S*“, the 


PA(w)(pg+( to)) is the available momentum of A(S+) 
in the center-of-mass frame of S*('0(3686)) at mass to, 
p(y(Ps+) is pK{rn)(p^+(rn)) for to = M, and L is the or¬ 
bital angular momentum. Due to the limited statistics, 
we do not determine the spin-parities of S(1690)“ and 
S(1820)“ with this data sample. In the fit, the spin- 
parities of S(1690)“ and S(1820)“ are assumed to be 
= 1/2“ and = 3/2“ based on previous experimen¬ 
tal results [7, 8 ], the angular momenta (^( 3 .-=;+)) 

are set to be 0 for both the S(1690)“ and S(1820)“, while 
the K~A angular momenta (T(a'-a)) ^^e 0 and 2 respec¬ 
tively. Bl{p) is the Blatt-Weisskopf form factor [22]: 


Bo{p) = 1; B2{p) = 



13 


3p2g2 


9gr 


( 3 ) 


where Qq is a hadron ’’scale” parameter which is on the 
order of 1 fm [22], and was set to be 0.253 GeV/c in the 
fit according to the result of the FOCUS experiment [23]. 

The overall fit result and the background components 
from the fit are shown as the solid and dashed curves in 
Fig. 3, respectively. The resulting masses, widths and 
event yields, as well as the corresponding significances 
of the S(1690)“ and S(1820)“ signals, are summarized 
in Table I, where the significance is evaluated by com¬ 
paring the difference of log-likelihood values with and 
without the S“ (1690/1820) included in the fit and tak¬ 
ing the change of the number of degrees of freedom into 
consideration. The significance is calculated when study¬ 
ing the systematic uncertainties sources (Sect. V) and 
the smallest value is reported here. The resonance pa¬ 
rameters from the PDG [2] are also listed in Table I for 
comparison. 

Due to the limited statistics, the measurement of spin- 
parity of S(1690/1820)“ is not performed in this anal¬ 
ysis. To determine the product branching fractions of 
the cascade decay 13(0(3686) —S(1690/1820)“S+) x 
13(S(1690/1820)“ — K~A), the corresponding detection 
efficiencies are evaluated with MC samples taking the 
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spin-parity of S(1690)“ and S(1820) to be = 1/2“ 
and Sj2~, respectively. The detection efficiencies and the 
corresponding product branching fractions are also listed 
in Table 1. Corresponding systematic uncertainties are 
evaluated in Sect. V. 

TABLE I. H(1690)” and H(1820)” fit results, where the 
first uncertainty is statistical and the second systematic. 
The B denotes the product branching fraction B{tp{3686) —> 
H(1690/1820)“H+) x Z3(H(1690/1820)“ R-A). 



H(1690)“ 

H(1820)“ 

M(MeV/c^) 

1687.7±3.8±1.0 

1826.7±5.5±1.6 

r(MeV) 

27.1±10.0±2.7 

54.4±15.7±4.2 

Event yields 

74.4±21.2 

136.2±33.4 

Significance)^) 

4.9 

6.2 

Efiiciency(%) 

32.8 

26.1 

B (10“®) 

5.21±1.48±0.57 

12.03±2.94±1.22 

MpDG(MeV/c^) 

1690±10 

1823±5 

rpDG(MeV) 

<30 

9^+15 

^^-10 


IV. ANALYSIS OF %l^{3686) ~fK-AE.+ 

In this analysis, the same selection criteria as those 
used in the 7/(3686) —>■ analysis are implemented 

to select the K~ and to reconstruct A and S"*" candidates. 
Photon candidates are reconstructed from isolated show¬ 
ers in EMC crystals, and the energy deposited in the 
nearby TOF counters is included to improve the photon 
reconstruction efficiency and the energy resolution. A 
good photon is required to have a minimum energy of 25 
MeV in the EMC barrel region (| cos(0)| < 0.8) and 50 
MeV in the end-cap region (0.86 < |cos(0)| < 0.92). A 
timing requirement (0 < t < 700 ns) is applied to further 
suppress electronic noise and energy deposition unrelated 
to the event. The photon candidate is also required to 
be isolated from all charged tracks by more than 10 °. 

The selected photons, K~, and A and candi¬ 
dates are subjected to a 4C-fit under the hypothesis of 
7/(3686) —>■ 7 Ar“AS+, and xic < 1^0 required. For 
events with more than one good photon, the one with 
the minimum xic i® selected. MC studies show that the 
background arising from 7/(3686) —>■ Ar“AS+ can be ef¬ 
fectively rejected by the 4C-fit and the xic requirement. 

With the above selection criteria, the M{ptt~) distribu¬ 
tion is shown in Fig. 4 (a). The A is observed clearly with 
low background, and the requirement 1.110 < M{p'!T~) < 
1.121 GeV/c^ is used to select A candidates. After that, 
the distribution of M(p 7 r+) is shown in Fig. 4 (b), where 
the A is observed with almost no background. The re¬ 
quirement 1.110 < M(p 7 r+) < 1.121 GeV/c^ is further 
applied to improve the purity. The M(A 7 r+) distribu¬ 
tion of the surviving events is shown in Fig. 4 (c), and 
a mass window requirement 1.315 < M(A 7 r+) < 1.330 
GeV/c^ is used to select 7/(3686) —?> 'yK~AE^ candi¬ 
dates. Figure 4 (d) shows the scatter plot of M{xA) 


versus M{K~A'EA) with all above selection criteria. 
The vertical band around the mass is from the de¬ 
cay 7/(3686) —>■ Ar“E°S+, while three horizontal bands 
around the XcJ (•^ = 0 , 1 , 2 ) mass regions are from 
7/(3686) lXcJ,XcJ -t K~AE+ (J = 0, 1, 2). There is 
also a horizontal band around the 7/(3686) mass region, 
which is background from 7/(3686) —>■ K~A'EA' with a 
random photon candidate. 

A. STUDY OF 7/(3686) ^ 

After applying all above selection criteria, the projec¬ 
tion of M{jA) is shown in Fig. 5, where a clear peak 
is visible with low backgrounds. As shown in Fig. 4 (d), 
the cascade process of 7/(3686) —>■ 7 Xc 2 , Xc 2 —t K~AE'^ 
will overlap with the band on M{xA). This pro¬ 
cess is investigated as potential background using the in¬ 
clusive MG sample together with the exclusive process 
7/(3686) —>■ 7 r+ 7 r“ J/t/, J/t/ —>• K~pYP. Both processes 
have the same final states as the signal, but do not pro¬ 
duce a peak in the M ( 7 A) distribution around the re¬ 
gion. The distribution of background obtained from the 
inclusive MG sample is shown as the shaded histogram 
in Fig. 5. The background is also studied with the can¬ 
didate events within the A or S"*" sideband regions of 
data, and the lack of peaking background in the M ( 7 A) 
distribution is confirmed. The background from e“'"e“ 
annihilation directly is estimated by imposing the same 
selection criteria on the data sample taken at y/s = 3.65 
GeV. No event survives, and this background is negligi¬ 
ble. 

To determine the 7/(3686) —>■ Ar“S°S+ yield, an ex¬ 
tended unbinned maximum likelihood fit of the M ( 7 A) 
distribution is performed with a double Gaussian func¬ 
tion for the together with a first order Chebychev 
polynomial for the background shape. The overall fit re¬ 
sult and the background component are shown in Fig. 5 
with solid and dashed lines, respectively. The fit yields 
the number of events to be 142.5±13.0, and the re¬ 
sulting branching fraction is ,6(7/(3686) —>■ A'“E°S+) = 
(3.67 ± 0.33) X 10“^, by taking the detection efficiency 
of 9.0% obtained from MG simulation and the branching 
fractions of intermediate states [2] in consideration. The 
errors are statistical only. 

B. STUDY OF XcJ IC“AH+ (J = 0, 1, 2) 

The XcJ —t A'“AS+ {J = 0, 1, 2) yields are de¬ 
termined by fitting the invariant mass distribution of 
K~ A'EA, M{K~ A’EA). To remove the background from 
7/(3686) — >■ K~YP’EA', the additional selection M{xA) > 
1.21 GeV/c^ is applied. The M{K~A'EA) distribution is 
shown in Fig. 6 , where the Xcj peaks are observed clearly. 
Potential backgrounds are studied using the events in 
the A or sideband regions of data and the inclusive 
MG samples. The inclusive MC M{K~AEA) distribu- 





a 

*► 


150 - 


100 - 



1.09 


1.10 


1.11 


1.12 


1.13 


1.14 


M(p 71 -) (GeV/cO 





FIG. 4. The invariant mass distributions of (a) p'K~, (b) (with the A selected) and (c) Att^. Dots with error bars are 
data, and the solid histogram is from the phase-space MC, which is normalized to the data. The arrows indicate the selection 
requirements used in the analysis (see text), (d) The scatter plot of M('yh) versus M{K~K'S^) for data. 



FIG. 5. The M{'yA) distribution, where the dots with error 
bars are data, the shaded histogram is the background con¬ 
tribution estimated from the inclusive MG sample, and the 
solid and dashed lines are the fit results for the overall and 
background components, respectively. 


tion is shown in Fig. 6 as the shaded histogram. Ac¬ 
cording to the MC study, the dominant backgrounds are 
from the cascade decays y/>(3686) —>■ J/ijj, J/tp —>• 

pK~TP —>■ yA and '0(3686) —^ , but 


none of them produce peak in the XcJ regions. 



FIG. 6. The K~AE'^ mass distribution, where the dots with 
error bars are data, the shaded histogram is the background 
contribution estimated from the inclusive MC sample, and 
the solid and dashed lines are the overall and background 
component contributions from the fit. 

An extended unbinned maximum likelihood fit of the 
M(Ar“AS+) distribution is performed to determine the 
number of XcJ events. The XcJ resonances are de- 
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scribed by Breit-Wigner functions convoluted with Gaus¬ 
sian functions to account for the mass resolution, and 
the background is described by a first order Chebychev 
polynomial function. The fit results are shown as the 
solid curve in Fig. 6, and the yields of Xcj {J = 0, 1, 2) 
are 56.9±8.9, 48.5±7.4 and 50.8±7.8 events, respectively. 
Taking the detection efficiencies, 6.9%, 8.5% and 6.9% 
for XcJ {J = 0, 1, 2) estimated by MC simulation, 
and the branching fractions of the decays of intermedi¬ 
ate states [2] into consideration, the product branching 
fractions B{^{3Q86) —>■ xXcj) x B{xcj iF“AS+) are 
measured to be (1.90 ± 0.30) x 10-^ (1.32 ± 0.20) x IQ-'^ 
and (1.68 ± 0.26) x 10“^ for Xcj {J = 0, 1, 2), respec¬ 
tively. The errors are statistical only. 


V. SYSTEMATIC UNCERTAINTY 

The different sources of systematic uncertainties for 
the measurement of branching fractions are considered 
and described below. 

a. Tracking efficiency 

In the analysis, both the proton and pion are from 
long lived particles (A or S“), and the corresponding 
tracking efficiencies are studied using a clean A control 
sample, selected by requiring the invariant mass recoiling 
against the pK^ system to be within the A mass region 
in the decay J/ip ^ ApK~^. The invariant mass recoiling 
against the pK^TT~ (pK^p) system is further required 
to be within the proton (pion) mass region to improve 
the purity of the control sample. The uncertainty of the 
tracking efficiency is estimated by the difference between 
efficiencies in data and MC samples and is parameter¬ 
ized as a function of transverse momentum. The average 
uncertainty of the proton (pion) tracking efficiency is es¬ 
timated to be 1% (1%) by weighting with the transverse 
momentum distribution of the signal. The uncertainty of 
the tracking efficiency is studied with a clean control 
sample of J/ip K*{892fK'^s, K*{892f and 

the systematic uncertainty is estimated to be 1% [24]. 

b. PID efficiency 

Similarly, the PID efhciencies of p/p and are esti¬ 
mated using the same control samples as those in tracking 
efficiency studies. All tracks are reconstructed and the 
target one is allowed to be unidentified. The systematic 
uncertainties for p, p and are all found to be 1%. 

c. Photon detection efficiency: 

The photon detection efficiency is studied utilizing the 
control samples ip{3Q8Q) —>■ tt^tt~J/iP, J/ip —>• p°7r° and 
ip{3&8&) —>■ Tr^TT^J/f/) with J/ip —>• {I = e, p) and 

The corresponding systematic uncertainty is esti¬ 
mated by the difference of detection efficiency between 
data and MC samples, and 1% is assigned for each pho¬ 
ton [25]. 

d. The secondary vertex fit: 

The efficiencies of the secondary vertex fits for A and 
are investigated by the control samples J/ip ^ ApK^ 
and J/ip ^ The differences of efficiencies between 


data and MC samples are found to be 1%, and are taken 
as the systematic uncertainties. 

e. Kinematic fit 

The track helix parameters {(po, k, tan A) for MC sam¬ 
ples are corrected to reduce the difference of the xic dis¬ 
tributions between data and MC [26]. The corresponding 
correction factors for kaons and the tracks from A de¬ 
cay (proton and pion) are obtained from a clean sample 
J/ip ^ ApK^, and those for the tracks from decay 
are obtained from the sample J/ip —>■ The sys¬ 

tematic uncertainties related to the 4C-fit, 1%, are esti¬ 
mated by the difference of efficiency between MC samples 
with and without the track helix parameter corrections. 

TABLE II. Summary of the relative systematic uncertain¬ 
ties (in %) in the branching fraction measurements. Here 
KAS,, KYPs,, XcJ and denote ip{3686) —>■ K~AS,'^, 

V>(3686) —>• A'“E°E+, ip{3686) 7Xcj, XcJ ^ K~A'EA and 
V>(3686) H(1690/1820)'H+,H(1690/1820)~ R-A, re¬ 

spectively. 


Source 

KAE 

KY'^E 

XcJ 


Tracking 

6 

6 

6 

4 

PID 

3 

3 

3 

3 

A vertex fit 

1 

1 

1 

1 

H vertex fit 

1 

1 

1 

- 

Kinematic fit 

1 

1 

1 

- 

Photon detection 

- 

1 

1 

- 

Signal model 

2.1 

0.5 

1.1,3.0,2.4 

0.8,1.6 

Background shape 

1.6 

0.5 

0,1.5,0.6 

7.1,7.! 

Fit range 

1.6 

1.9 

0.2,0.1,0.2 

6.3,4.4 

Mass shift, resolution 

- 

- 

- 

0.6,0.4 

Mass windows 

2.9 

1.4 

3.2,2.3,1.8 

1.0,1.3 

B{A — >• p-TT”) 

0.8 

0.8 

0.8 

0.8 

B{E+ Atv+) 

0.035 

0.035 

0.035 

0.035 

B{A —>• p7r+) 

0.8 

0.8 

0.8 

0.8 

B{ip{3686) 7Xcj) 

- 

- 

3.2,4.3,4.0 

- 

A^7(3686) 

0.8 

0.8 

0.8 

0.8 


Total 8.2 7.6 8.5,9.3,8.7 11.0,10.1 


/. The fit method: 

The systematic uncertainties related to the fit method 
are considered according to the following aspects. 
(1) The signal line-shapes: In the measurements of 
B{ip{3686) —>■ A'“AS+), B{xcj —>■ Ar“AS+) and 
B{ip{3686) —>■ Ar“S°S+), the signal line-shapes are re¬ 
placed by alternative fits using MC shapes, and the 
changes of yields are assigned as the systematic uncer¬ 
tainties. In the measurements of B{ip{3686) —>■ S*“S+), 
the corresponding uncertainties mainly come from the 
uncertainty of Qq. Alternative fits varying the Qo values 
within one standard deviation [23] are performed, and 
the changes of yields are treated as the systematic uncer¬ 
tainties. (2) The background line-shapes: In the measure¬ 
ments of B{ip{3686) —>■ Ar“AS+), B{xcj —>■ AT"AS+) and 
B{ip{3686) —>■ A'“E°S+), the background shapes are de¬ 
scribed with a first order Chebychev polynomial function 
in the fit. Alternative fits with a second order Cheby¬ 
chev polynomial function are performed, and the result¬ 
ing differences of the yields are taken as the systematic 
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TABLE III. Summary of the systematic uncertainties on H* parameters. 



H(1690) 
M (MeV^2) 

r (MeV) 

H(1820) 
M (MeV/c^) 

r (MeV) 

Signal model 

0.2 

0.3 

1.5 

1.2 

Background shape 

0.3 

1.8 

0.5 

3.3 

Fit range 

0.3 

1.7 

0.2 

2.2 

Mass shift, resolution 

0.5 

0.8 

0.2 

0.2 

Mass windows 

0.7 

0.7 

0.4 

0.6 

Total 

1.0 

2.7 

1.6 

4.2 


uncertainties related to the background line-shapes. In 
the measurement of B{ip{3686) —>■ an alterna¬ 

tive fit with a reversed ARGUS function (rARGUS) *, 
FrARGVsirn) = UARGUs(2mo - to), for the non-resonant 
components is performed, where toq is the mass thresh¬ 
old of K~A. The changes in the yields are taken as sys¬ 
tematic uncertainties. (3) Fit range: Fits with varied fit 
ranges, i.e., by expanding/contracting the range by 10 
MeV/c^ and shifting left and right by 10 MeV/c^, are 
performed. The resulting largest differences are treated 
as the systematic uncertainties. (4) Mass shift and reso¬ 
lution difference: In the measurement of branching frac¬ 
tions related to a Gaussian function G{m, 

which represents the mass resolution, is included in 
the fit, where the parameters of Gaussian function are 
evaluated from MC simulation. To estimate the system¬ 
atic uncertainty related to the mass shift and resolution 
difference between data and MC simulation, a fit with a 
new Gaussian function with additional parameters, i.e., 
G(m, fj, -I- Aii,a -I- Act), is performed, and the resulting 
difference is taken as the systematic uncertainty. The 
additional values A/i and Act are estimated by the dif¬ 
ference in the fit results of the A and between data 
and MC simulation. 

g. Mass window reguirement: 

The systematic uncertainties related to A and mass 
window requirements are estimated by varying the size 
of the mass window, i.e. contracting/expanding by 2 
MeV/c^. The resulting differences of branching fractions 
are treated as the systematic uncertainties. 

h. Other: 

The systematic uncertainties of the branching fractions 
of the decays ■i/'(3686) —>■ 7Xcj, —>• Att” and A —>• pTT~ 

are taken from the world average values [2]. The uncer¬ 
tainty in the number of f(S686) events is 0.8%, which is 
obtained by studying inclusive ip{3686) decays [11]. The 
uncertainty in the trigger efficiency is found to be negli¬ 
gible due to the large number of charged tracks [27]. 

The different sources of systematic uncertainties in the 
measured branching fractions are summarized in Table II. 
Assuming all of the uncertainties are independent, the 


* The ARGUS function is defined as TargusI”^; ™o, c,p) = m(l — 
-expfcfl — where mo is the mass threshold and 

c and p are parameters fixing the shape 


total systematic uncertainties are obtained by adding the 
individual uncertainties in quadrature. 

In the measurement of the resonance parameters, 
the sources of systematic uncertainty related to the fit 
method and the A and S'*" mass window requirements are 
considered. The same methods as those used above are 
implemented, and the differences of the mass and width 
of are regarded as the systematic uncertainties and 
are summarized in Table III. The total systematic uncer¬ 
tainties on resonance parameters obtained by adding 
the individual uncertainties in quadrature are shown in 
Table III. 


VI. CONCLUSION 

Using a sample of 1.06 x 10® tp{3686) events collected 
with the BESIII detector, the processes of tp{3686) —>■ 
K~ KEF and tp{3686) — >■ ^K~ MEA are studied for the 
first time. In the decay i/'(3686) —>■ A'“AS+, the branch¬ 
ing fraction ,B('0(3686) —>■ K~KE'^) is measured, and 
two structures, around 1690 and 1820 MeV/c^, are ob¬ 
served in the K~ A invariant mass spectrum with signif¬ 
icances of 4.9 ct and 6.2 ct, respectively. The fitted reso¬ 
nance parameters are consistent with those of S“(1690) 
and S“(1820) in the PDG [2] within one standard devi¬ 
ation. The measured masses, widths, and product decay 
branching fractions B{tp{3686) — >■ S*“S+) x B{E*~ —>• 
K~A) are summarized in Table I. This is the first 
time that S~(1690) and S~(1820) hyperons have been 
observed in charmonium decays. In the study of the 
decay -0(3686) —>■ ^K~AE^, the branching fractions 
B{tp{3686) —>■ Ar“S°S+) and B{xcj —t K~AE'^) are 
measured. All of the measured branching fractions are 
summarized in Table IV. The measurements provide new 
information on charmonium decays to hyperons and on 
the resonance parameters of the hyperons, and may help 
in the understanding of the charmonium decay mecha¬ 
nism. 
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TABLE IV. Summary of the branching fractions measurements, where the first uncertainty is statistical and the second 
systematic. 


Branching fraction 


Decay 


(3.86 ±0.27 ±0.32) x 10"'’ 
R-K (5.21 ± 1.48 ±0.57) X 10"® 
JfA (12.03 ±2.94 ±1.22) x 10"‘ 

/o I r\ oo I r\ r»r>\ 


i/)(3686) K-AS.+ 
V>(3686) H(1690)-H+, =(1690)” 

V>(3686) ^ H(1820)“H+, H(1820)~ 
V>(3686) ^ K-Ty=.+ 

i>{S686) jxco, Xco R-AS."' 

^(3686) 7Xci, Xci R~AE'^ 

^(3686) 7Xc2, Xc2_-> R~AE"' 

Xco ^ K-AE+ 

Xci ^ R-AE+ 

Xc2 ^ R-AE+ 


(3.67 ±0.33 ±0.28) x 10”® 
(1.90 ±0.30 ±0.16) X 10"® 
(1.32 ±0.20 ±0.12) X 10"® 
(1.68 ±0.26 ±0.15) X 10"® 
(1.96 ±0.31 ±0.16) X 10"^ 
(1.43 ±0.22 ±0.12) X lO""' 
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